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Abstract 

We derive an asymptotic solution of the vacuum Einstein equations 
that describes the propagation and diffraction of a localized, large- 
amplitude, rapidly-varying gravitational wave. We compare and con- 
trast the resulting theory of strongly nonlinear geometrical optics for 
the Einstein equations with nonlinear geometrical optics theories for 
variational wave equations. 

1 Introduction 

Geometrical optics^ and its generalizations, such as the geometrical theory 
of diffraction, are a powerful approach to the study of wave propagation, 
for both linear and nonlinear waves. In this paper, we develop a theory 
of strongly nonlinear geometrical optics for gravitational wave solutions of 
the vacuum Einstein equations. Specifically, we derive asymptotic equa- 
tions that describe the diffraction of large-amplitude, rapidly-varying grav- 
itational waves. 



^Here, we use the term 'geometrical optics' to refer to any asymptotic theory for the 
propagation of short- wavelength, high-frequency waves, irrespective of its area of physical 
application. 
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These equations are a generalization of the straightforward non-diffractive, 
nonhnear geometrical optics equations for large-amplitude gravitational waves 
derived in [2|. This strongly nonlinear theory differs fundamentally from the 
weakly nonlinear theory for small-amplitude gravitational waves obtained by 
Choquet-Bruhat and Isaacson jT^, because it captures the direct nonlin- 
ear self- interaction of the waves. 

The Einstein equations may be derived from a variational principle, and, 
when written with respect to a suitable gauge, they form a system of wave 
equations for the gravitational field. In order to explain the structure of 
nonlinear geometrical optics theories for gravitational waves, and to moti- 
vate the form of our asymptotic expansion, it is useful to consider first such 
theories for a general class of variational wave equations. 

We describe straightforward and diffractive geometrical optics theories 
for variational wave equations in Section |2 In the weakly nonlinear the- 
ory for waves with periodic waveforms, the amplitude- waveform function 
a{6, v) of the wave depends on a ray variable v and periodically on a 'fast' 
phase variable 6. The wave-amplitude a satisfies the Hunter-Saxton equa- 
tion (EH, 

|a, + (^^«')^ + ^«}^ = l^i'^l - («e>} , (1-1) 

where the angular brackets denote an average with respect to 9. The coef- 
ficient A of the nonlinear terms in may be interpreted as a derivative 
of the wave speed with respect to the wave amplitude. In addition, the 
wave-amplitude is coupled with a slowly-varying mean field. 

An application of this expansion to the Einstein equations leads to the 
theory of Choquet-Bruhat |3] and Isaacson In that case, the nonlinear 
coefficients corresponding to A in the equations for the amplitudes of the 
gravitational waves are identically zero. This theory therefore describes 
a nonlinear interaction between a high-frequency, oscillatory gravitational 
wave and a slowly-varying mean gravitational field, but does not describe 
the direct nonlinear self-interaction of the gravitational wave itself. 

In Section [2n we make a distinction between 'genuinely nonlinear' wave- 
fields in variational wave equations, for which A is never zero, and 'linearly 
degenerate' wave-fields for which A identically zero (see Definition^. All 
wave-fields in the Einstein equations are linearly degenerate. As observed in 
[HI [Zl ini 1 for example, this fact reflects a fundamental degeneracy in 

the nonlinearity of the Einstein equations in comparison with general wave 
equations. 
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A property of the Einstein equations related to their hnear degeneracy 
is that they possess an exact solution for non-distorting, large-amplitude, 
plane waves, the Brinkman solution |H1 119j . One can use this solution to 
derive a strongly nonlinear geometrical optics theory for large-amplitude 
gravitational waves. We outline the resulting non-diffractive theory [2 in 
Section E31 

We derive our diffractive, strongly nonlinear geometrical optics solution 
of the vacuum Einstein equations in Section ^ This solution describes the 
propagation and diffraction of a thin, large-amplitude gravitational wave, 
such as a pulse or 'sandwich' wave. The simplest, and basic, case is that 
of a plane-polarized gravitational wave diffracting in a single direction. We 
summarize the resulting asymptotic equations here. 

We suppose that the polarization of the wave is aligned with the diffrac- 
tion direction. Then, with respect to a suitable coordinate system (u, v, y, z), 
the metric g of the wave adopts the form 



g = -2e-^^^{du - eYdy)dv + e'^ie^dy^ + e'^dz'^) + 0{e'^). (1.2) 



Here, e is a small parameter. The leading-order metric component functions 
{U, V, M), and the first-order function Y depend upon a 'fast' phase variable 
6, an 'intermediate' transverse variable rj, and the 'slow' ray variable v, where 



The phase u is light-like and the transverse coordinate y is space-like. 

This metric describes a gravitational wave whose wavefronts are close 
to the null-hypersurface u = 0. The wave is plane-polarized in the (y, z)- 
directions and diffracts in the y-direction. 

To write equations for the metric component functions in a concise form, 
we define a derivative Df^ and functions (j), tJj by 



where (| denotes the 'raising' operator from one-forms to vector fields. 




(1.3) 
(1.4) 
(1.5) 



DrjM-e'^Ye 
D„{U + V). 



We note that D^i has the geometrical significance that 



I?, = -e^(dy)» 
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Then {U, V, M, Y) satisfy the following system of PDEs: 



Uee - \ {Ui + Vg^) + UeMe = 0, (1.6) 
{cl) + iP)e = m + V)e, (1.7) 
Ue, - UeU, = ie~(^+^+^^) |z)^</. + - ^(^^ - H - ^fj^^ (1.8) 

Ve, - \ {UeV, + U,Ve) = ie-(^+^+^^) | + ^0'} , (1-9) 
Me. + ^ ([/,[/. - VeV.) = le'^^+v+M) |_^^^ _ 1^2 ^ ^aj ^ (^^^g) 

This system is the main result of our analysis. It is an asymptotic reduction 
of the full vacuum Einstein equations to a (1 + 2)-dimensional system of 
PDEs. The system provides a model nonlinear wave equation for general 
relativity, and should be useful, for example, in studying the focusing of 
nonplanar gravitational waves and the effect of diffraction on the formation 
of singularities. We plan to study these topics in future work. 

We also derive asymptotic equations for the diffraction of a gravitational 
wave of general polarization in two transverse directions. These equations 
are given by (|3lT|) - (|3T3|) . (jHH) . and (jn^ - lfriKI) . but they are much more 
complicated than the ones for plane-polarized waves written out above. 

We study the structure of (|1.6|) - (|l.lflj) in SectionlH The constraint equa- 
tion (|1.6() is a nonlinear ODE in 6 in which 77, v appear as parameters. As 
shown in Proposition this constraint is preserved by the remaining equa- 
tions. Equation H1.7() is a linear, nonhomogeneous ODE in 9 for y, whose 
coefficients depend on 9- and ry-derivatives of (C/, M). It may therefore 
be regarded as determining Y in terms of (J7, V, M). Equations (|1.8jl - (|l.l()|) 
are a system of evolution equations for (?7, V, M) which are coupled with 
Y . The main part of the system of evolution equations consists of a (1-1-2)- 
dimensional wave equation in {9,r],v) for (U + V), in which {9,v) are char- 
acteristic coordinates. 

When Y = and all functions are independent of rj, equations ()1.6|1 - 
(|1.10|) reduce to the colliding plane wave equations, without the v-constraint 
equation (see Section . When all functions depend on (C,r/) with ^ = 
9 — Xv for some constant A, we get a system of PDEs in two variables that is 
studied further in [T] . This system describes space-times that are stationary 
with respect to an observer moving close to the speed of light. 
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2 Geometrical optics 



The Einstein equations do not form a hyperbolic system of PDEs because 
of their gauge-covariance, but their properties are closely related to those of 
hyperbolic PDEs. In order to develop and interpret geometrical optics solu- 
tions for the Einstein equations, it is useful to begin by studying geometrical 
optics solutions for hyperbolic systems of variational wave equations. 

We remark that there is a close correspondence between geometrical 
optics theories for hyperbolic systems of conservation laws |12| and varia- 
tional wave equations. For example, the inviscid Burgers equation [31 El 
is the analog of the Hunter-Saxton equation (|2.14|) . and the unsteady tran- 
sonic small disturbance equation jllj is the analog two-dimensional Hunter- 
Saxton equation (|2.15|1 . One can also derive large-amplitude geometrical 
optics theories for linearly degenerate waves in hyperbolic conservation laws 
(see (231 121] ! for example) that are analogous to the large-amplitude theories 
described here for the Einstein equations. 

There are other nonlinear geometrical optics theories for dispersive waves, 
most notably Whitham's 'averaged Lagrangian method' |28l I29j for large- 
amplitude dispersive waves. This theory has a different character from the 
ones for nondispersive hyperbolic waves. Nonlinear dispersive waves have 
specific waveforms — given by traveling wave solutions — in which the 
effects of dispersion and nonlinearity balance, whereas nondispersive hy- 
perbolic wave equations, and the Einstein equations, have traveling wave 
solutions with arbitrary waveforms, which may distort as the waves propa- 
gate. 

In straightforward theories of geometrical optics, waves are locally ap- 
proximated by plane waves and propagate along rays. In linear geometrical 
optics, the wave amplitude satisfies an ODE (the transport equation) along 
a ray; in nonlinear geometrical optics, the amplitude- waveform function typ- 
ically satisfies a nonlinear PDE in one space dimension (a generalization of 
the transport equation) along a ray. This difference may be understood as 
follows: linear hyperbolic waves propagate without distortion, so one re- 
quires only an ODE along each ray to determine the change in the wave 
amplitude; by contrast, wave-steepening and other effects typically distort 
the waveform of nonlinear hyperbolic waves, so one requires a PDE along 
each ray to determine the change in the wave amplitude and the waveform. 

Straightforward geometrical optics, and a local plane- wave approxima- 
tion, break down when the effects of wave-diffraction become important; 
for example, this occurs when a high-frequency wave focuses at a caustic, or 
when a wave beam of large (relative to its wavelength), but finite, transverse 
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extent spreads out. The effects of diffraction on a high-frequency wave may 
be described by the inclusion of additional length-scales in the straightfor- 
ward geometrical optics asymptotic solution. 

Perhaps the most basic asymptotic solution that incorporates the ef- 
fect of wave diffraction leads to the 'parabolic approximation', described 
below for wave equations. In Section |31 we will derive analogous asymptotic 
solutions of the Einstein equations for the diffraction of large-amplitude 
gravitational waves. 



2.1 The wave equation 

We begin by recalling geometrical optics theories for the linear wave equation 

9tt = y- [clVg) . (2.1) 

Here, g{t,x) is a scalar function, the wave speed CQ{t,x) is a given smooth 
function, and x G R*^. Although well-known, these theories provide a useful 
background to our analysis of nonlinear variational wave equations and the 
Einstein equations. 

We look for a short-wavelength asymptotic solution g = oi (|2.1|) . 
depending on a small parameter e, of the form 

, , , / uit, x) \ 
5f(t, x)~al 2 — j as e ^ 0. 

The solution depends on a 'fast' phase variable^ 

and the 'slow' space-time variables (t,x). 

One finds that the scalar-valued phase function n(t, x) satisfies the eikonal 
equation 

= cglVtip. (2.2) 
The amplitude-waveform function a{6,t,x) satisfies the equation 

ae^ + Nae = 0, (2.3) 

where 

= Utdt - clVu ■ V (2.4) 



■^Here, we use u/e'^ as a phase variable, rather than u/e, for consistency with the 
difFractive expansion below. 



6 



is a derivative along the rays associated with the phase u, and N{t, x) is 
given by 

N = ^{uu-V-{clVu)}. (2.5) 

Equation H2.3|l has solutions of the form 

a{e,t,x) = A{t,x)F{e) , 

where A{t, x) is a wave-amplitude, which satisfies an ODE along a ray (the 
transport equation of linear geometrical optics) 

A^ + NA = 0, (2.6) 

and F{6) is an arbitrary function that describes the waveform of the wave. 
For example, if F{9) = e*^, then the solution describes an oscillatory har- 
monic wave; if 

r e>o, 

6* < 0, 



F{9) 



then the solution describes a wavefront across which the normal derivative 
of g of order n jumps; and if F[9) has compact support, then the solution 
describes a localized pulse. 

The term A'^^ in the transport equation (|2.6|) describes the effect of the 
ray geometry on the wave amplitude. The coefficient becomes infinite 
on caustics where the rays focus. The straightforward geometrical optics 
solution breaks down when this happens, and diffractive effects must then 
be taken into account. 

There are many ways in which diffraction modifies straightforward geo- 
metrical optics. Here, we consider one of the simplest diffractive expansions, 
given by 

^u {t,x) y{t,x ) 

9 ' 



g^{t,x)'-^a[ 2 — ) ,t,xj as e ^ 0. 



This asymptotic solution depends upon an additional 'intermediate' variable 



y 



where y{t, x) is a scalar-valued transverse phase. 

One finds that u satisfies the eikonal equation, as before, and 

Vv = 0, 
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meaning that y is constant along the rays associated with u. Moreover, the 
ampUtude-waveform function a(^, 77, t, x) satisfies the equation 

aev + Nae + ]^Da^^ = 0, (2.7) 

where the coefficient D{t, x) of the diffractive term is given by 

D = yl-cl\Vy\\ (2.8) 
For harmonic solutions, we have 

a{9,r],t,x) = A{r],t,x)e^^, 
and equation 1)2.71) reduces to a Schrodinger equation 

i{A^ + NA} + ^DA^^ = 0. 

This 'parabolic approximation' and its generalizations are widely used in 
the study of wave propagation [27] . 

In the simplest case of the diffraction of plane wave solutions of the 
two-dimensional wave equation with wave speed cq = 1, 

gtt = 9xx + Oyy , 

we may choose 



u = —^, y = y, 



t + x 



V2' 

Equation (|2.7|) is then 

This equation describes waves that propagate in directions close to the x- 
direction, and is a wave equation in which 6 and v are characteristic coordi- 
nates. We will see that equations with a similar structure to ()2.9p in their 
highest-order derivatives arise from the Einstein equations. 

2.2 A variational wave equation 

Next, we consider the following nonlinear, scalar wave equation Jl] 

fftt - V • {c\g)Vg) + c{g)c'{g)\Vg\^ = 0. (2.10) 
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We assume that the wave speed c : M ^ M'*' is a smooth, non-vanishing 
function, and a prime denotes the derivative with respect to g. This equation 
is derived from the variational principle 

^ I {lst-lc\9)\Vg\^^ dtdx = 0. 

The structure of the nonlinear terms in (|2.in|) resembles that of the 
Einstein equations, although, as we shall see, the effects of nonlinearity 
are qualitatively different because of a 'linear degeneracy' in the Einstein 
equations. 

We look for an asymptotic solution of 1)2. 1U() of the form ^ ^1 
g (t, X) ~ go{t, X) + e a ( — ^ — ,t, x j as e — > 0. 

This solution represents a small-amplitude, high-frequency perturbation of 
a slowly- varying field go. The amplitude and the wavelength of the pertur- 
bation are chosen to be of the same order of magnitude because this leads 
to a balance between the effects of weak nonlinearity and the ray geometry. 

First, we suppose that the amplitude-waveform function a{9, t, x) is a 
periodic function of the phase variable 0. We assume, without loss of gen- 
erality, that its mean with respect to 9 is zero. 

We find that the phase u{t,x) satisfies the linearized eikonal equation 
(|2.2|1 . with Co = c{go). The mean- field go satisfies the nonlinear wave equa- 
tion 

gott - V • {clVgo) + coc'o\Vgof + \Vufcoc'o{al) = 0, (2.11) 

where the angular brackets denote an average with respect to 9 over a period, 
and Cg = c' (go)- Equation (|2.11|) has the same form as the original wave 
equation, with an additional source term proportional to the mean energy 
of the wave- field. 

The amplitude- waveform function a{9, t, x) satisfies the periodic Hunter- 
Saxton equation (|1.1() . 

|a, + (^Aa2^^ + iVaj^ = ^Aja^ - {aj}} , (2.12) 

where dy is the ray derivative defined in 1)2. 4() . N is given by 1)2. 5|) . and 

A = -IVuI^coCq. (2.13) 

The coefficient A is proportional to the derivative of the wave-speed with 
respect to the wave amplitude, and it provides crucial information about 
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the effect of nonlinearity on the waves; for a given go, we obtain a nonhnear 
PDE for the amphtude-waveform function only if A 7^ 0. 

The expansion summarized above for waves with periodic waveforms is 
uniformly valid for {t,x) = 0(1) and u = 0{e~^) as e — > 0. If we consider 
localized waves, then we may obtain an asymptotic expansion that is valid 
near^ the wave-front n = 0, in which we neglect the mean-field effects. Thus, 
the background field goix) is a solution of the original wave equation, and 
a{6, t, x) satisfies 

|a„+QAa2^ +iVa| = ^Aal (2.14) 

It is straightforward to include diffractive effects in the expansion for 
localized waves. The asymptotic solution has the form 

, , , 9 / u(t,x) v(t,x) \ 
g^t,x) ~ go{x)+e^a ( ^ ^ as e ^ 0. 

Here, go is a solution of the original wave equation, the phase u satisfies the 
eikonal equation at go, the transverse function y is constant along the rays 
associated with u, and a{6,r],t,x) satisfies a (1 -|- 2)-dimensional generaliza- 
tion of the Hunter-Saxton equation 1)2. 14(1 . 

|a„ + Q^a^^ + Naj + ^Da^^ = ^Aaj, (2.15) 

where D is given by (|2.8jl . 



2.3 Systems of variational wave equations 

In this section, we consider a class of hyperbolic systems of nonlinear wave 
equations that are derived from variational principles of the form Jl] 

^j^v'^9)^^,dx = ^. (2.16) 

Here, x = {x^ , . . . , x*^) G R'^^^ are the space-time variables, 

g= {g\...,g'^) : R^+l ^ 

are the dependent variables, Apq : R™ — > R are smooth coefficient functions, 
and we use the summation convention. We assume that Apg = Apq = Aqp . 

Specifically, the expansion is valid when 6 — 0{1), or u = 0{e^), and t = 0(1). 
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The Euler-Lagrange equations associated with ()2.16|1 are 



Gpb]=0, (2.17) 

where 

We assume that ()2.17() forms a hyperbohc system of PDEs. The scalar wave 
equation considered in the previous section is the simplest representative of 
this class of equations. 



The weakly nonlinear geometrical optics solution of (|2.17j) has the form 

PHI EI 

g^'ix) goix) + e^a(^^^,x^ R{x) as e ^ 0, (2.19) 

where go : M'^^^ M"* is a slowly-varying function, u : M'^+i ^ M is a 
phase function, o : M X M'^+i ^ R is an amplitude-waveform function, and 
R : M'^+^ — > is a suitable vector field. We find that the phase u satisfies 
the eikonal equation 

det C = 0, (2.20) 
where the m x m matrix C{x) has components 

Cpg = u,<^u,pA^^{go). (2.21) 

The vector R= {R^,..., R'^Y in (|2l9|) is a right null- vector of C, so that 
CpqR'^ = 0. Here, and in ()2.19|) . we assume that we are dealing with a simple 
characteristic, meaning that the null-space of C is one-dimensional. 

We write the scalar-valued amplitude-waveform function a{6,x) as a 
function of the 'fast' phase variable 9 = u/e"^, and the slow variables x. 
If a{9, x) is a periodic function of 9 with zero mean, then the mean-field 
goix) satisfies the equation 

Gp[go] = lHp{al), (2.22) 
where the angular brackets denote an average with respect to 9, and 

a/3 



8 A' 

Hp = u^.u,,^{go)Rmr 
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Equation H2.22() has the same form as the original equation for g with an 
additional source term proportional to (flg). 
The equation for a is 1)2. 12() with 



5„ = 2u^,A'^^{go)Wmd^c,, 

A = u,<^u,,^{go)RPRm\ (2.23) 

For localized solutions, the mean- field interactions may be neglected, and 
we can obtain diffractive versions of this expansion as before. 

These equations generalize easily to the case when equation 1)2. 17(1 has 
a multiple characteristic of constant multiplicity n > 2, say. In that case, 
one obtains a mean-field equation whose source term is a sum of averages 
of products of ^-derivatives of the wave amplitudes, and an n x n coupled 
system of Hunter-Saxton equations for the wave amplitudes. In particular, 
if the vectors {Ri, ■ ■ ■ , Rn} form a basis of the null-space of the matrix C 
defined in (|2.21|) . then the coefficients Aijk of the nonlinear terms in the 
system of equations for the wave amplitudes are given by 

^ijk = Ux'^u,,-g^{go) R^iR'Rl l<i,j,k< n. (2.24) 

2.4 Linearly degenerate wave equations 

Motivated by the corresponding definition for hyperbolic conservation laws 
introduced by Lax ^01, we make the following definition. 

Definition 1 A simple characteristic of a hyperbolic system of wave equa- 
tions i2.17\ )- i2.1f^) is genuinely nonlinear (respectively, linearly degenerate^ 
if, for every go G M™ and every non-zero du G M'^+^ that satisfies \2.2(l\) . 
the quantity A defined in \2.2'J\) is non-zero (respectively, zero). A multiple 
characteristic of constant multiplicity n is linearly degenerate if all coeffi- 
cients Aijk defined in \2.24^ are zero. We say that the system is genuinely 
nonlinear ( respectively, linearly degenerate ) if all of its characteristics are 
genuinely nonlinear (respectively, linearly degenerate). 

Thus, for linearly degenerate wave equations, the amplitude- waveform 
function in the weakly nonlinear theory corresponding to the ansatz (|2.19)) 
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always satisfies a linear PDE. Interactions with a mean-field may still occur, 
however. 

For example, from H2.13() . the nonlinear wave equation ()2.1U() is genuinely 
nonlinear if c'{g) ^ for all g €M, and linearly degenerate if c is constant, 
when it reduces to the linear wave equation. 

Even if a wave- field is not linearly degenerate, a loss of genuine nonlin- 
earity may occur at a particular point G M."^ and direction du G M*^^"^. 
A loss of genuine nonlinearity is related to the null condition for nonlinear 
wave equations introduced by Klainerman jl7j . The Einstein equations do 
not satisfy the Klainerman null condition completely because of the non- 
zero coupling between gravitational waves and a mean gravitational field 
[S]. Nevertheless, it is the vanishing of the nonlinear A-coefficients in the 
equations for the amplitude-waveform functions at the Minkowski metric, 
for example, that permits the existence of global smooth, small-amplitude 
perturbations [HI El l22]- A similar result would not be true for (1 + 3)- 
dimensional variational wave equations in which A 7^ 0. 

We emphasize that this definition for wave equations is analogous to but 
different from the definition for hyperbolic systems of conservation laws. If a 
genuinely nonlinear wave equation were rewritten as a first-order hyperbolic 
system, it would be classified as a linearly degenerate first-order hyperbolic 
system, since its wave speeds are independent of the derivative Vg. 

Hyperbolic systems of conservation laws with linearly degenerate wave- 
fields possess non-distorting plane wave solutions, and one can derive a large- 
amplitude geometrical optics theory for them |231 124j . Linearly degenerate 
variational systems of wave equations possess non-distorting plane wave so- 
lutions, and a large- amplitude geometrical optics theory, only if they satisfy 
certain additional degeneracy conditions. We will not analyze these issues 
here, however, and instead discuss the analogous theory for the Einstein 
equations. 

2.5 The Einstein equations 

The Einstein equations do not fall exactly into the class of variational wave 
equations considered in Section [2.31 because of their gauge-covariance. They 
can be derived from a variational principle of the form 1)2. 16() . obtained after 
an integration by parts in the Einstein-Hilbert action, but the resulting 
Euler-Lagrange equations are not hyperbolic. Nevertheless, as is well-known, 
they become hyperbolic when written with respect to a suitable gauge, such 
as the harmonic gauge. The equation for the metric g, with components g^^ 
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and determinant g, then adopts a similar form to (|2.17|) - H2.18|1 . namely 



where Hf^i, is a quadratic form in the metric derivatives dg with coefficients 
depending on g. (See for example, for an explicit expression.) 

The weakly nonlinear expansion described in Section [2.31 for variational 
systems of wave equations corresponds to the expansion of Choquet-Bruhat 
[1] and Isaacson jl5j for the Einstein equations. The Einstein equations are 
not strictly hyperbolic, and all of the nonlinear coefficients in the asymptotic 
equations for the amplitude-waveform functions are zero. Thus, in this 
generalized sense, the Einstein equations are a linearly degenerate system of 
variational wave equations. 

The Einstein equations possess a non-distorting, plane wave solution for 
large-amplitude gravitational waves, which forms the basis of a strongly 
nonlinear geometrical optics theory for large-amplitude gravitational waves. 
It does not, however, appear possible to obtain a self-consistent theory for 
oscillatory large- amplitude waves, since the mean energy- momentum associ- 
ated with an extended wave-packet would generate a very strong background 
curvature of space-time. This restriction is related to the fact that mean- 
field interactions with oscillatory waves already occur at leading order in the 
small-amplitude theory. We therefore consider localized waves. 

An asymptotic theory for the propagation of localized, large-amplitude, 
rapidly varying gravitational waves into slowly varying space-times was de- 
veloped in PJ. In this theory, the metric of a plane-polarized wave may be 
written, with respect to a suitable coordinate system {u,v,y,z), as 



Here, the metric component functions {U,V,M), depend on {0,v), where 
the 'fast' phase variable 6 is given by 




The functions {U, V, M) satisfy the following PDEs: 




g = -2e-*^dn dv + e'^ie^dy^ + e'^dz^) + O(e^). 



(2.25) 




(2.26) 



- UgU, = 0, 

- \ {UeVy + U,Ve) = 0, 
+ I {UeU, - VeV,) = 0. 




'V 



(2.28) 



(2.29) 
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Equations (|2.27|) - H2.29|) are wave equations for (U, V, M) in characteristic 
coordinates {0,v), and H2.26() is a constraint which is preserved by ()2.27|) - 

These equations correspond to a well-known exact solution of the Ein- 
stein equations, the colliding plane wave solution ^1^1123120]) without 
the usual constraint equation in v, 

Uvv - \ ipl + K') + U^M, = 0. 

This constraint need not be satisfied if the slowly-varying space-time into 
which the wave propagates is not that of a counter-propagating gravitational 
wave. If it is not satisfied, then the resulting metric is an asymptotic solution 
of the Einstein equations but not an exact solution. 

The asymptotic equations we derive in this paper are a diffractive gen- 
eralization of (IT^ - (fT^ . 



3 The asymptotic expansion 

In this section, we outline our asymptotic expansion of the Einstein equa- 
tions, and specialize it to the case of plane-polarized gravitational waves 
that diffract in a single direction. Our goal is to explain the structure of the 
expansion and the resulting perturbation equations. The detailed algebra is 
summarized in the appendices. 



3.1 The general expansion 

Let g be a Lorentzian metric. We denote the covariant components of g 
with respect to a local coordinate system by ga/3- The connection coeffi- 
cients T^a/s and the covariant components Rafs of the Ricci curvature tensor 
associated with g are defined by 



2 \dx°' dxf^ ^x^' , 
" ^~d^ ~ ^ r^a/jT^AM - ^'"axT^iSt^- (3.2) 

Here, and below, Greek indices a,f3,X,fi run over the values 0,1,2,3. The 
vacuum Einstein equations may be written as 

Rap = 0. (3.3) 
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We look for asymptotic solutions of ()3.3() with metrics of the form 



g = g f x"; ej , (3.4) 

gie, r/^ e) =g (9, rf , x") + e g (0, rf , x") + g (0, r,«, x") + ©(e^). 

Here, e is a small parameter, u is a phase, y° with a = 2, 3 are spacelike 
transverse variables, and 

^ "w r, y"" 
= -, v^ = — 

e 

are 'stretched' variables. The ansatz in 1)3. 4jl corresponds to a metric that 
varies rapidly and strongly in the w-direction, with less rapid variations in 
the y'^-directions, and slow variations in x". 

We remark that the form of this ansatz is relative to a class of coordinate 
systems, since an e-dependent change of coordinates can alter the way in 
which the metric depends on e. It would be desirable to give a geometrically 
intrinsic characterization of such an ansatz, and to carry out the expansion 
in a coordinate-invariant way, but we do not attempt to do so here. 

Using ()3.4() in the Einstein equations ()3.3() . we find that, to get a non- 
trivial solution at the leading orders, the phase u must be null to leading 
order, and the transverse variables must be constant along the rays as- 
sociated with u to leading order. Moreover, the leading-order metric must 
have the form of the colliding plane wave metric. 

To carry out the expansion in detail, we use the gauge-covariance of the 
Einstein equations to make a choice of coordinates that is adapted to the 
metric in (|3.4|) . We assume that u is approximately null up to the order 
e^, that y'^, are constant along the rays associated with u to the order e. 
Furthermore, we assume that we can extend (n, y'^,y^) to a local coordinate 
system 

(x°,x^x2,x^) = {u,v,y'^,y^). 

Then, as shown in Appendix El we can use appropriate gauge transfor- 
mations, which involve a near identity transformation of the phase and the 
transverse variables, to write a general metric (|3.4j) whose leading-order term 
has the form of the colliding plane wave metric, as 

g = 2 5-01 dx"dx'^+ Qab dx^dx^ 

+e |2 Qia dx^dx'^+ 9ab dx^dx^^ + hj dx'dx^ + 0{e^). (3.5) 
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In H3.5() . and below, indices a,b,c, . . . take on the values 2, 3, while indices 
i,j,k,... take on the values 1,2,3. We raise and lower indices using the 
leading order metric components; for example, we write 



a/3 





9 aP 



We have one additional gauge freedom in (|3.5() . involving a nonlinear 
transformation of the phase, which we can use to set either 





9oi= 1 



or one of the three components 



112 

9i2, 9n, 9ii 



equal to zero. We will exploit this gauge freedom later, since it is conve- 
nient in formulating a variational principle for the asymptotic equations (see 
Appendix Id]). 

Using the method of multiple scales, we expand derivatives of a function 
/ with respect to x'^ as 



£2' e 



e 



e 



(3.6) 



and then treat (0, ry"^, x) as independent variables. In (|3.6|) . in the appendices, 
and below, we use the shorthand notation 



f =^ 



77" ,x 



du 
dxi^ 



_ dy°- 
dxi^ 



df 



dxi^ 



We use (|3.5|) - (|3.6|1 in (|3.H1 - (|3.2|1 and expand the result with respect to 
e. We find that 



r a/3 + - r a/3+ r af3 + 0{e), 

1-4 1-3 1-2 _^ 

^"/3 ~ ^ +^ RafS +^ i?a/3 +0{e ) 



(3.7) 



where explicit expressions for the terms in the above expansions are given 
in Appendix^ The nonzero components of the Ricci tensor up to the order 
e~2 are summarized in Appendix O 
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Using 1)3 .711 in ()3.3() and equating coefficients of e ^, e ^ and e ^ to zero, 
we get 

0, (3.8) 

0, (3.9) 

0. (3.10) 

These perturbation equations lead to a closed set of equations for the leading- 
order and first-order components of the metric, as we now explain. 
The only component of 1)3. 8|) that is not identically satisfied is 

Roo= 0. (3.11) 
From (|3.11|) and (jC'.ll) . it follows that 9^/3 satisfies 

"2(5 9ab,e),e + 2 5 9oi,e9''^ 9ab,e 9^' 9bc,e9'"^ 9ad,e= 0. 

This equation is a constraint on the leading order metric which involves only 
^-derivatives, and has the same form as the constraint equation for a single 
gravitational plane wave. 

The nonzero components of the Ricci tensor at the next order in e are 

-3 -3 

Roo, Roa ■ 

The condition 

Roa= 0, (3.12) 

1 

yields equations for 5ia and their derivatives with respect to 6, with non- 
homogeneous terms depending on 9- and r/-derivatives of the leading order 
metric. Thus, (|3.12)) provides two equations relating the ffi'st-order pertur- 
bation in the metric to the leading order metric. 
The equation 

i?oo=0 

1 

is a single equation that is homogeneous in the first-order components dab, 
as can be seen from ()(].2|) . A nonzero solution of this equation corresponds 
physically to a free, small-amplitude gravitational wave of strength of the 
order s propagating in the space-time of the large-amplitude gravitational 



-4 

RafS = 
-3 

Raf3 = 
-2 

RaB = 
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wave. Retaining a nonzero solution of this homogenous equation would not 
change our final equations for the large-amplitude wave, and for simplicity 
we assume that 

Sab= 0. 

This assumption is also consistent with the higher-order perturbation equa- 
tions. 

At the next order in e, the nonzero components of the Ricci curvature 

are 

-2 -2 -2 -2 

ROO, ROI, ROa, Rab ■ 

-2 -2 

The corresponding perturbation equations for Rqq and i?oa in (!3.10j> give 
non-homogenous equations for the second-order metric components 

2 2 
9ab, 9 la, 

with source terms depending on the lower-order components of the metric. 
These equations are satisfied by a suitable choice of the second-order metric 
components, and they are decoupled from the equations for the leading- 
order and first-order metric components. We therefore do not consider them 
further here. 

The remaining perturbation equations are 

~Roi= 0, 1^ab= 0. (3.13) 

By use of the gauge freedom mentioned at the beginning of this section, we 
may set 

2 

9ii= 0. 

In that case, (|3.1H|) provides a set of four equations relating the leading-order 
and first-order metric components, and their derivatives. 

We remark that, in general, one cannot eliminate secular terms from the 
asymptotic solution that are unbounded as the 'fast' phase variable 9 tends 
to infinity. As a result, the validity of the asymptotic equations is restricted 
to a thin layer of thickness of the order about the hypersurface u = 0, 
where 9 = 0{1). A global asymptotic solution can be obtained by matching 
the 'inner' solution inside this layer with appropriate 'outer' solutions, such 
as slowly varying space-times on either side of the wave. In this paper, 
however, we focus on the construction of asymptotic solutions for localized 
gravitational waves, and do not consider any matching problems. 
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Summarizing these results, we find that the asymptotic solution is given 

by 

g = 2 5oi dx^dx^+ 9ab dx^dx^ + e |2 5ia dx^dx"^ + 0{£^), (3.14) 
where the six metric components 

1 

9m, 9ab, 9la 

satisfy a system of seven equations (|3.1ip - ()3.13p . Equation ()3.11|) is an ODE 
in 0, and, in the case of plane-polarized waves, we show that it is a gauge- 
type constraint that is preserved by the remaining equations H3.12() - ()3.13() . 

The explicit form of the equations follows from the expressions in ()C.1|) , 
(|(].3|) ~ (|(].5|) for the corresponding components in the expansion of of the 
Ricci tensor that appear in H3.11() - (|3.13() . In general, these equations are 
very complicated, but they can be simplified considerably in special cases. 



3.2 Plane-polarized gravitational waves 

In this section, we specialize our asymptotic solution to the case of a plane- 
polarized gravitational wave that diffracts in a single direction, and write 
out the resulting equations explicitly. 
We choose coordinates 

x\ x^) = (u, V, y, z) (3.15) 

in which the metric has the form (|3.5j) . We suppose that the metric depends 
on the variables {9,r],v), where the phase variable 6 and the transverse 
variable rj are defined by 

and is independent of the second transverse variable C = z/e. This means 
that the wave diffracts only in the y-direction. We could also allow the metric 
to depend on z, which would appear in the final equations as a parameter. 

We consider a leading-order metric that has the form of the collid- 
ing plane wave metric for a plane-polarized wave, polarized in the (y, z)- 
directions. We have seen that we need only retain the higher order compo- 
nents 

1 2 
9la, 9ii . 



20 



One component of H3.12() is homogeneous in 5i3, since there is no dependence 
on C, and without loss of generahty we can take 

9i3= 0. 

We therefore take the special form of the metric (|H.5)) given by 

g = -2e-*^ (^du - eYdy - ^^^Tdv^ dv + e'^ie^dy^ + e'^dz'^) + 0{e'^), 

(3.16) 

where the functions {U,V, M,Y,T) depend on (6,ri,v). We recall that we 
have the gauge-freedom to set either M, Y, or T equal to zero. In writing 
the equations, we choose to set T = 0. 

The metric 1)3. 16() must satisfy equations ()3.11() - (|3.13|) . First, using 
(|3.16|) in (|3.11|) . and simplifying the result, we obtain, after some algebra, 
the ^-constraint equation ()1.6() . 

Next, using (|3.16|1 and ()(].3|) . we find that the only nontrivial component 
of (pn^ is 

Ro2= 0. (3.17) 

Furthermore, after the introduction of (j), ip defined in (|1.4|) - (|1.5|) and some 
algebra, we find that (|3.17|) may be written as equation (|1.7() . 

Finally, using (|3.1(ij) and (|(].4|) ~ (|(].5|) . we find that the only nontrivial 
components of 1)3. 13(1 are 

iloi=0, R22=0, R33=0- (3.18) 

After some algebra, these equations may be written as (jl.8|) ~ (jl.in|) . where 
the derivative is defined in H1.3() 

3.3 The variational principle 

Equations (j3.1H) - (j3.13() can be derived from a variational principle, which is 
obtained by expanding the variational principle for the Einstein equations. 
In particular, the variational principle provides a check on the algebra used 
to derive the equations. (We have also checked the final results by use of 
MAPLE.) 

In order to formulate a variational principle, it is necessary to retain the 
component 

2 

9ii . 
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Variations with respect to this component yield the gauge-type constraint 
H3.11() . and it may be set to zero after taking variations with respect to it. 

The general form of the asymptotic variational principle is given in Ap- 
pendix^ Here, we specialize it to the case of a plane-polarized gravitational 
wave considered in Section 13.21 Using the metric H3.16|l in (|D.4() we find, 
after some algebra, that the variational principle HD.5|) for ()1.6|) - H1.1U|) is 

SS^^^ = 0, 5(1) = j dOdvdr], 

where the Lagrangian L^"^) may be written as 

L^-^^ = e-^{2Me, + AUg, - VgV, - 3UeU, 

-Tee + Te{Me + 2Ue) + T{Mee + 2Uee - pi - \vi) 

^^-iU+V+M) [_2^^^ _ ^ 3^2 ^ 2</>V; + ^^]]. (3.19) 

Variations of 5^^) with respect to the second-order metric component T lead 
to the ^-constraint equation (|1.6I) . Variations with respect to the first-order 
metric component Y , which appears in tp, and D,,, lead to equation (|1.7() . 
Variations with respect to f7, V ^ M lead to the evolution equations ()1.8|1 - 
(frrnil . after we set T = 0. 

4 Properties of the equations 

In this section, we study some properties of the equations for plane-polarized 
waves that we have derived in the previous section. We write out the struc- 
ture of the highest-order derivatives that appear in the equations, and use 
this to formulate a reasonable IBVP for them. We also show that the 6- 
constraint equation is preserved by the evolution in v, and that the linearized 
equations are consistent with the linearized equations for gravitational waves 
in the parabolic approximation. 

4.1 Structure of equations and an IBVP 

In this subsection, we consider the structure of equations (|1.6j) - (|l.lfl|) in 
more detail. The first equation, ()1.6() . is an ODE with respect to 6 relating 
{U,V,M). As we show in the next section, this equation is a gauge-type 
constraint which holds for all v if it holds for u = 0, say. We may therefore 
neglect this equation provided that the initial data at u = is compatible 
with it. 
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The remaining equations (|1.7|1 - H1.10|) form a system of equations for 
(U,V, M,Y). In order to exhibit their structure, we rewrite them in a way 
that shows exphcitly how the highest, second-order, derivatives appear. 

Using ()1.3() - (|1.5() . we may rewrite equation H1.7|) as 

Yee - {{V + M)eY}g + |^[/| - UgVe - WgA Y 



= {U + V + M)g^ + M^Ue -{U + V\Ve. (4.1) 

If (f7, V, M) are assumed known, then this equation is a hnear ODE in for 
Y. 

We may rewrite equations (|1.8I) - H1.10|) as 

{U + V + M)e.-\{U + V),{U + V)^ 

= ie-(^+^+^^)|-l02-,^^}, (4.2) 
{U + V)e, - \ue (U + V), - hj, {U + V), , 

(4.4) 

Examining the terms that involve second order derivatives, we see that these 
equations consist of a (1 + 2)-dimensional wave equation in {9, rj, v) for (U + 
V), and two (1 + l)-dimensionaI wave equations in (9, v), for (U + V + M) 
and U, in which {9, v) are characteristic coordinates. An additional second- 
order derivative term, proportional to Dj^Yg, appears in the equation for U. 
The function Y is also coupled with the evolution equations through the 
dependence of the transverse derivative Dr/, given in (|1.3j) . on Y. 

In order to specify a unique solution of these equations, we expect that 
we need to supplement the ODE (|4.1|) with data for Y and Yg on 9 = 0, 
say, and the evolution equations (|4.2() - (|4.4() with characteristic initial data 
for {U, V,M) on9 = and ^; = 0. Thus, a reasonable IBVP for (|IT|) - (|01l 
in the region 9 > 0, — oo < r] < oo, and v > is 

{U,V,M) = {Uo,Vo,Mo) onv = 0, 
{U,V,M) = {Ui,Vi,Mi) on 9 = 0, 
iY,Yg) = {Yo,Yi) on 9 = 0. 
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Here, {Uo,Vo, Mq) are given functions of {0,r]) that satisfy the constraint 

Uoee - \ {U^e + Voe) + UoeMoe = 0, 

and (Ui,Vi, Mi,Yq,Yi) are given functions of {r],v). This data may be 
interpreted as initial data for the state of the wave on the hypersurface 
f = 0, and boundary data on the leading wavefront = 0. 

4.2 The constraint equations 

In this section, we show that the constraint equation (|1.6j) is preserved. 

Proposition 1 Suppose that (U, V, M, Y) are smooth functions that satisfy 
[T7^- in~m} . Let 

F = Ue9-\ {Ui + Ve^) + UgMg. (4.5) 

Then 

= U,F. (4.6) 



Proof. We write the evolution equations H1.8|) - H1.1U|1 as 

Uo. - UeU. = ie-(^+^+^^)^, (4.7) 

Ve, - \ {UeV, + U,Ve) = ie'^^+^+^^^S, (4.8) 

Me, + \ {UeU, - V,Ve) = I^-^u+v+m)^^ ^^ g^ 

where, using the notation defined in H1.3() - ()1.5|l . 

^ = D^</. + D^V-^'^'-'/'V'-^', (4.10) 

B = -Dr,(t> + \<l}\ (4.11) 

C = -D^^l,-}-ct? + ilj^. (4.12) 

Differentiating 1)4. 5|) with respect to f , and using H4.7|) - (|4.9|) and ()4.5p to 
replace Uq,, Vgv,Mg, and Ugg in the result, we find that 

= C/,F + ie-(^+^+^^) A (4.13) 
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where 

D = Ae-{U + V)gA-VeB + UeC. (4.14) 

Differentiating equation (|4.1UI) for A with respect to 9, introducing the com- 
mutator [dg,Dr^] of do and D^, and using equations ()1.5|) and (|1.7() . we 
compute that 

Ae = deS^Dr^i^ + i;)-^^^ -4>i;-i;^^ 

= Dr,{ct> + i^)e + [de, Dr,] + V) - + i^)4>9 -{4> + 2V)^9 
= n,{i^{U + V)e} 

-{4> + ij)(f>e -(</> + 2V)Ve + [de, Dr,] (</> + tp) 
= {U + V)eDr,i; + tPDr,{U + V)e 

-{4> + i;)(f>e -{((> + 2V')Ve + [de, D^] + ^) 
= iU + V)eDr,il> + i^deD^{U + V)-^ [de, {U + 

-(</> + -(</> + 2V)V'e + i^r,] (<A + V') 
= ([/ + F)eI?^V' + i^i^e -(</> + - (0 + 2V')^0 

+ [ae, (0 + V') - ^ [a^, D.,] (u + y) 

= {U + V)eDr,i;-{cl) + i;){ct> + i')e 

+ [de, n,] (0 + V) - [de, D,,] {U + V) 

= {U + V)e{Dr,i;-i;{(l) + i:)} 

+ [de, n,] + V-) - [de, Dr,] {U + V). (4.15) 

From the definition of D,, in (|1.3|) . it follows that 

[de,Dr^] = UeDr^ + e^Yede. 

Hence, using H1.4() - ()1.5() and H1.7|) . we get 

[de,Dr,]{^ + ^P) = UeDr,{(l) + i>) + e^Ye{^ + ^)e 

= UeDr,{(P + iP)+i;e^Ye{U + V)e, 

[de,Dr,]iU + V) = UeDr,{U + V) + e^Ye{U + V)e 

= Ue^ + e^YeiU + V)e. 

Using these equations in (|4.15j) . and simplifying the result, we find that 

Ae = iU + V)e {D^,^ - cP^; - V^^) + Ue {D^^cP + Dr,^; - V'') ■ 

Finally, using this equation and ()4.1U() - ()4.12() in 114. 14(1 . and simplifying the 
result, we find that D = 0. It follows from that F satisfies □ 
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4.3 Linearization 

We consider the small-amplitude limit of (|1.6jl - p.lUj) in which 

U, V,M,Y ^ 0. 

From the constraint equation ()1.6|) . we have U = 0{V'^), so U is of higher 
order in a linearized approximation and can be neglected completely. From 
H1.3|) - H1.4|) . we also have in this approximation that 

D — 8 

<j) = M^-Ye, 

^ = 

Linearization of H1.7|) yields 

+ = 0, 

while linearization the evolution equations H1.8() - (|1.1U() gives 

^ev = ~2'^V' 
Mev = -^'^v 

Neglecting some arbitrary functions of integration for simplicity, we find 
that these equations are satisfied if M = —V and Y = 0. In that case, V 
satisfies an equation of the form (|2.9j) : 

1 

The corresponding linearized metric is given by 

g = 2(1 - V)dudv + (1 + V)dy^ + (1 - V)dz'^. 

Linearization of the Einstein equations, with a suitable choice of gauge, 
leads to a set of linear wave equations for the metric components. One can 
verify that this linearization of our asymptotic solution agrees with what 
is obtained by an application of the parabolic approximation, described in 
Section 12.11 to the linearized Einstein equations. 
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A Gauge transformations 



We consider a Lorentzian metric g, and local coordinates x", in which 

g = Oa^dx" dxf^. 

We denote the contravariant form of g by g". 

We suppose that the metric g depends upon a small parameter e, and 
there exist independent functions y"(x), a = 2,3, such that as e — > 

we have 

g«(dn,(fu) = 0(e2), (A.l) 
^^du^dy") = 0{£), a = 2, 3. (A.2) 

The first condition states that du is an approximate null form up to the 
order e^. The second condition states that is approximately constant 
along the rays associated with u. 

In this appendix, we reduce such a metric to the standard asymptotic 
form, given in (|3.5)) . It is convenient to consider two different types of 
coordinate transformations. First, we show that it is possible to make a 
near-identity change of coordinates, to a local coordinate system x° with 

such that 

n 1 h 

g = 2 5oi dx"dx'^+ gab dx^-dx" 

+£ |2 gia dx^dx"'+ lab + % dx'dx^ + O(e^). (A.3) 

Second, we show that the form of the metric can be further simplified by 
use of a coordinate x^ such that 

xO = u + 0{e'^). 

As before, indices a, f3, . . . take on the values 0, 1, 2, 3; indices a,b,c, . . . take 
on the values 2, 3; and indices k, . . . take on the values 1, 2, 3. 
Let x" be any local coordinate system in which 

x° = u, x'' = y". 

Then the conditions (|A.1|) . HA.2|) are equivalent to 

aOo4°° = 0, 9^- = 0. (A.4) 
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These conditions imply immediately that 

9u= 0, 5ia= 0. (A.5) 

Moreover, we have 

det{gab) > 0, (A.6) 

which says that are space-like independent coordinates. 

The expansion (|l-{.4j) . with the additional restriction ()A.5|) . is invariant 
under a coordinate transformation of the form 



x^ + e^^-"! ^,^,:E;e) , (A.7) 



^\9, r?^ x; e) = e ^0 (0, r,", x) + |,o (e, i^\x) + ©(e^), 

^^(0, r/^ x; e) =1' (6, r?", x) + e (0, r/^ x) + (0, x) + 0{e^), 



where, as usual, 



^0 

— , = — 



To begin with, we perform a change of coordinates of the form 

x°^x°, ^ x' + ^\e,7]^,x). (Ai 

At order zero in e, the metric components transform according to 

foo ^ 5oo+2*|'0 5ofc + ^'|'e^',e5«, 

0,0 



9ij ^ 

Recalling HA.6|1 . the matrix (gab) is invertible, and we can choose 

. 

so that 

9oo= 0, 9oa= 0. (A.c 

^ ^ 

Thus, the zero order metric g takes the form given in ()A.3|) . 
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Next, we consider the effect of the coordinate transformation HA.7|) on 
the metric components at order e. Under the action of the transformation 

^ + h"{9,r]^,x), (A.IO) 

the form ()A.3|) of the metric is unchanged at order zero. At order one, the 
components transforms according to 

1 1 1 -I 

900 5oo+2^,e9oi, 

1 1 1 „ 

901 9oi + ^"ggoi, 

1 1 1,0 

90a 9oa + '^,e9ab, 

1 1 

9ij ^ 9^,. 

1 

Choosing appropriately ^f"^, these transformations can be used to make 

90a= 0. (A.ll) 
Moreover, the condition (|A.4() imphes 

in=-(5oi) 5°° = 0. (A.12) 

Next, we simphfy the metric components at order e^. Under the action 
of the transformation 

^ + |"(0,r/^x), (A.13) 

the form of the metric is unchanged at order zero and one. At order two, 
the metric components transform according to 

2 2 2 

900 foo+2^^0foi, 

2 2 2 

901 5oi + ^°0 5oi, 

2 2 2 

90a 5oa + ^|'e5a6, 

2 2 
9^3 ^ 

2 

Choosing x]/" appropriately, these transformations can be used to make 

9oa= 0, (A.14) 
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and the metric can be written as in ()A.3|) . 

If we allow for changes of order in x", then we still have the freedom 
to use a transformation of the form 

^ e^-^(Q,rf,x). (A.15) 

The substantial difference between HA.13|) and HA.15() is that, unlike the for- 
mer one, the latter involves a nonlinear transformation of the phase Q. This 
transformation can be interpreted as a transformation to characteristic co- 
ordinates in the space (0, ry'*, x). The gauge ()A.15|) acts only on the following 
components: 

oil 022 

ffoi^ 5oi, 9ia^9ia +^',a foi, 9u^9u +2^,e 9oi . 

We can use ()A.15|) to set ^01= 1 or one of the three components 

112 
9i2, 5i3, 9ii 

equal to zero. More generally, it is possible to choose ^ so that an appro- 
priate relationship involving all of these components is satisfied. 
The contravariant form of the metric tensor in (|A.3|) is 

g« = 2 5°iaoai+ g'^'dadb - e |2 9'^''doda+ g^'dadb 

+ 2 h^''doda+ h^'dadbj + 0(8^), (A.16) 
with 





_(|oo_ 




lody 




2n 


^Oc 


lad^ 






1 





In l|A.16|) . 5°^ and 9"-^ are defined by 

Oni 00 

5°^ 9oi= 1, 5"^^ 9cb= 51 

o„, , 

Moreover, 9 is used to raise the sub-index 1, and 9 is used to raise the 

sub-indices 2,3. Thus, for example, we have 

l00^(001)2 1^^^ loajg.l^gabl^^^ 

It follows from HA.16|) that the transformed phase u = x^ and transverse 
variables y°- = x°- satisfy ()A.H) - (|A.2|) . 
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B Expansion of connection coefficients and Ricci 
tensor 



We look for a multiple-scale expansion as e — > of the metric components 
of the form 

/ u{x) y"'{x) \ 

gap = 5q/3 I ^2~'^^'^!^ ) ' ^'^■^1 

(e, 7/^ x; e) = {9,r]'' ,x) + elap {9,^ ,x) + l^p {9 .rj" ,x) + 0{e^). 
The contravariant metric components satisfy 

Expansion of this equation in a power series in e gives 

:= 5"^-e5"^-e2 (5"^- V^^^S'^M +0(e3), (B.2) 

„ . . 

where 9 is the inverse of 9a(3- In (|B.2|) and below, we use the leading 

order metric components to raise indices, so that 

g^P =g"^^ g^- g^^, 9^^ =g"^' 9^" 9f,u ■ (B.3) 

With this notation, note that, for instance, the first order term in the ex- 
pansion of the covariant metric components (/"^ with respect to e is not 9"^ 

but = - ff"^. 

The expansion of the connection coefficients (or Christoffel symbols) is 



A 1 ^^A 1 ^^A ° A 

^ af3 = T a/3 + -r Q-/3+ F o/3 + 0{e), 



where 



-2, 1 Oa,; /^o 

r a/3 = -7: 9 [ 9i3fM,e Ua+ 9ati,e up- 9al3,e Uf, 



-1. I . / 

r\/3 = -9^^'[9p^,-aya+9af,,ayp-9aP,ay^ 



1 Oa„ /^l 1 1 

+ - 5' 5/3/,,£) Ua+ 9atJ.,e up- 9ai3,e u^ 
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-- 5 {9f3f,,e Ua+ 9at,,e Up- 9af3fi I , 

+ ^9^'' (9pf,,a ya+ 9ai,,a Vp- 9af3,a 

-- ^-^'^ i9f3^-a yl+ 9 at, -a Vp- 9af3,a < J 

1 , /2 2 2 \ 

+ - 5 I 5/3/,,0 Ua+ 9ap.,e Up- 9ap,e U^ I 

1 1, /l 1 1 \ 

-- 5 ( 5,3/i,e ^*Q+ 5Q^,e 9a/3,e 1 

1 2 . /O \ 

+- /i ( 9p^,e Ua+ 9ati,e up- g^pfi u^ \ . 
The expansion of the Ricci tensor is 

1-4 1-3 1-2 _-|^ 

^oifi = ^ Ra(3 +^ Rafi +^ i?Q/3 +0{e ), 

where 

-4 -2 -2 -2 -2 -2 -2 

Ra(3 = r ^a(3,eUf,— T ^PfM,eUa+ F '^a/S F '^^i/— T ^oi/ T 



V 

-3 -1 -1 -2 -2 

Rafi = r r ^Pp,fiUa+ r ''a/3,a?/;j- T ^/J/^.a^a 

—2 ~1 —1 —2 —2 —1 —1 —2 

i cxfi i /^i/T^ i ^ fiu L au 1- p^ i ai/ i 



-2 -1 -1 

Rap = T^aPfiU^,- T^Pf,fiUa+ T ^aP,ayl- T ^ Pf,,ayl 

-2 -2 -2 -2 



-2 -2 -2 -2 

-2 -2 -1 -1 -1 -1 

pM pJ^^-pM pJ^^-LpM^p!^ _pM pi 



C Nonzero connection coefficients and Ricci ten- 
sor components 

It this appendix, we write out expressions for the nonzero connection co- 
efficients and Ricci tensor components for a metric with the form given in 
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Nonzero connection coefHcients at order e ^: 

-2q Oqj^ -2-^ 1 Oqj^ -2^ 1 0^^ 

r 00 =9 9oi,0, r ab = "2 ^ ^aM' ^ 06 = 2 ^ ^6c,6l • 

Nonzero connection coefficients at order £~^: 

i Oa — 2 ^ (yOl.a + yia,e) ~ 2 9ab,0, 

r la=2^ {9oi,a - 9la,e), T ab = -29 9ab,d, 

r "01 = "2 ^"^(ffoi.c - ^ice), r "ob = - 5"'^ ^bce "2 ^"'^ ^foc,^, 

-1„ 10,0 , 1 1„ 

r "be = 2 9''\9m,-c + - + 2 56c,e • 

Nonzero connection coefficients at order e^: 

r 00 — y 501,0, r 01 — 2 ^ \9m,c — 9ic,e) + 2 ^ii,^'' 

On 10 m 2 llnl 12„0 

r%„ = - 9°\9la,e + 5oi,a) - 2 5ac,e +2 ^ 9ac,e, 

Oq 1 Oq^ 1 1 

r ab- 2^ Wl6,a + 9ia,b " ^a&.l j 

1 l0c/° ° ° ^ 1 ^00 ° 

, 1 0„-, 2^0, 1 0„, 2 ^ 

T\i = - 9^H2 <7oi,i - 9n,0), T\a = ^ 9°H9oi,a - 9ia,e), 

T^ab = - ^ 9^^{9ab,0 + 9ab,e), 

102 111 

r%i = - g^^'iOice - 9oi,c) - 2 ff"'(5ic,e - 5oi,c), 

102 111 120 

r"o6 = - 9"''^{9bc,e + 9bcfi) - 2 ^'"^ 9bc,e, +2 ^^<^fi' 

loo 1 1 llnO 1 

T\b = 2 9''''{9bc,i + 9,,;, - 9ib,c) - 2 5 " ^i^), 

llfil 12„0 10,1 1 1 

T%c = 2 9 9bc,e "2 9bc,e +2 5"'^(56d,c- + 5,^,6 - ^bcd) 

1 1 .0 0,10,0 , 
-2 5"''(56d,c- + 9,a,-b - 9 bed) + 2 5 (^Hc + 5cd,6 - <?6c,d)- 

Nonzero component of tfie Ricci curvature tensor at order e"^: 

-4 10,0 lOniO 0,0 lO 0,. 

Roo= "2(5 9ab,e),e + ^ 3 9oi,e9'''' 9ab,e 9"" 9bc,e9^'^ 9ad,e ■ (C.l) 
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Nonzero components of the Ricci tensor at order £ ^: 

-3 11 10 n, 1 10,, 1 

i?oo= -2 K,ee + 2 ^ a^ifidle - ^ 9'"^ 9bc,e9le, (C.2) 

-3 1 .0 Om 1(, . 10,0 Om 1;, 1 .0. 

Roa= 2(fa65 f + ^ 9"^ 9cd,e9ab9'^^ 9% + -(5"" 5a6,e),c 

1 .Oni 0,0 , 10,0 0.0 

-^(9 9oi,a + 9^ 9cd-a),e + 4 9ab,e9'^^ 9de,c 

1 Oni 0.0 1 0. , 0^„ 

+^ 5°' 9oi-a9'"' 9cd,e -4 9^^ 9cd,e9''' 9be,-a ■ (C.3) 

The nonzero components of the Ricci curvature at the order s are Rqo, 
-2 -2 -2 -2 -2 -1 

Roi, Roa and i^a^- The components Rqo and i?oa are the analogs of Roa 
-1 

and Roa for the higher order components of the metric, and will not be 
listed here. The remaining components of the Ricci tensor are listed below: 

-2 /OniO 10, \ 10,0 0,,0 

Roi= - (^5"' <7oi,i +2 5 9cd,iJ ^ - 4 ^ 9bc,e9'"' 9ad,i 

+ - l^i^Ol.a - 9la,e) 9 J ~ 2 y-901,a - 9la,0) 9 
10 1 /inO.O 0,0.0 \ 

+^{9oi,-a - 9ia,e) 9"' 9cd,e - 5"' 9"' ' (C-4) 

-2 1 0„, 0,0 ^ 

Rab= - 9^^ 9ab,W +2 9'"^{9ac,e9bd,l + 9ac,l9bd,e) 

10„, 0,0 * 

9^^ 9'^{9cd,l9ab,e + 9cd,e9ab,l)+ Rah 

1 ,ln , 1 In Om 1 , 

+ 2(5 9abfi)-c +^9^"^ 9ab,e (5°^ 901,c +-9'^^ 9de,c) 

+ \(9^^ 9ab,c),e + \ 9^" 9ab,c (5°^ 9oifi +\ 9'^" 9de,e) 
+-(5' 9ac9^p,0 + -9 9ac9^i[9 9Qifl+-9 9de,e) 

l.OmOl^, lOniOl^OniO lO.O 

9bc9'i,a),e + 2 ^ 9bc9la{9^^ 9oi,e +2 5 9de,e) 

1 Ool Ocd 1 1 Ocrf ^Oe .0 ^ n 

-- 5°^ ^^-^^ 5ac,e<?ife,J + 2 5 9ac,e {9be,d " ^ferf.e) 

--9''^9"' 9bc,99ia,d +2 5^ 56c,e " 9ad,e) 

l,lnlnO , InlnO .OniO lO.O 

--{9'i 9ab,e),e- 5° 9^' 9ab,e 5oi,e +2 5* 5de,e) 
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llnlnO.O 10 mO 1 0^0 1. 

+- 5° 5°^ 9'"' 9ac,e9M,e 9ac9le 5°' 5^^^?,^. (C.5) 

* 

Here, are the components of the Ricci tensor of the zero order metric 



9af3 regarded as a function of the variables rj^ only, 
1 fo^a x\ 

Rab- 2 {9bd,a + 9ad,b " 5'afe,dJ J _ 

+ 1 9"^{9bd,a + - 9oi,c +^ 9'^ 9ef,c) 

-{9 9oi,a +- 9 9cd,a)fi ~ 2 90l,a9 9oi,b 

1 0^. 

9'^{9de,a + 9,,^a - 9ad,e) 9''^ {9bf,c + 9,f-^ - 9,^j). 

D Variational principle 

The variational principle for the vacuum Einstein field equations is 

55 = 0, ^ = / 

L = R^J-detg, (D.l) 
where R is the scalar curvature, 

R = 9''^Rap. 

Using (|B.1() . HB.2() . and (|3.7|) to expand the scalar curvature, we obtain 

1-3 1-2 , 
R+^ R +0(e" ), 



R 


1 


-4 

+; 


-4 

R 




-4 


-3 

R 


J.. 


-3 


-2 

R 




-2 



(D.2) 



1 o -3 2^-4 

For a metric of the form ()A.3|) . (|3.1(ij) . with 9ii= T, we find that 

-4 -3 

R = R= 0, 

-2 Oni -2 0,-2 In -3 2 „„ -4 

R= 2 i^oi + 9"^ Rab -2 5°° iJoa + i?oo • (D.3) 
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We use (|D.3|) in ()D.1|1 and expand the result with respect to e. This 
gives 

L = lL(-2) + 0(e-i), 



£2^ 



with 



L(-2) = )2 Roi+ goig"'' Rab -2 5? iJoa + /i? i?oo [ V det gab- (D.4) 

We make the change of variables in the integration 

d^x = du dv dy dz = e^d9 dv dr] dz, 

and omit the integration with respect to the parametric variable z. The 
leading order asymptotic variational principle then becomes 

55^=0, S^^^ = J L^-^Uedvdr]. (D.5) 

Variations of S^^^ with respect to the second order metric component ^ii 
give the constraint H3.11() . Variations with respect to the first order metric 

1 

components ^oa give the equation H3.12|) . Variations with respect to ^oi and 



gab give the evolution equations H3.13|) . 
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